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Abstract 

The equations providing to describe generation process in FEL with varied parame- 
ters of diffraction grating (photonic crystal) are obtained. It is shown that applying 
diffraction gratings (photonic crystal) with the variable period one can significantly 
increase radiation output. It is mentioned that diffraction gratings (photonic crystal) 
can be used for creation of the dynamical wiggler with variable period in the sys- 
tem. This makes possible to develop double-cascaded FEL with variable parameters 
changing, which efficiency can be significantly higher that of conventional system. 
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1 Introduction 



Generators using radiation from an electron beam in a periodic slow-wave 
circuit (travelling wave tubes, backward wave oscillators, free electron lasers) 
are now widespread [1]. 

Diffraction radiation [2] in periodical structures is in the basis of operation 
of travelling wave tubes (TWT) [3,4], backward wave oscillators (BWO) and 
such devices as Smith-Purcell lasers [5,6,7] and volume FELs using two- or 
three-dimensional distributed feedback [8,9,10,11]. 

Analysis shows that during operation of such devices electrons lose their en- 
ergy for radiation, therefore, the electron beam slows down and gets out of 
synchronism with the radiating wave. These limits the efficiency of generator, 
which usually does not exceed ~ 10%. 

In the first years after creation of travelling wave tube it was demonstrated 
[4] that to retain synchronism between the electron beam and electromag- 
netic wave in a TWT change of the wave phase velocity should be provided. 
Application of systems with variable parameters in microwave devices allows 
significant increase of efficiency of such devices [4,12]. 

The same methods for efficiency increase are widely used for undulator FELs 
[13]. 

In the present paper we consider generation process in Smith-Purcell FELs, 
volume FELs, travelling wave tubes and backward wave oscillators using pho- 
tonic crystal built from metal threads [14,15,16,17]. It is shown that applying 
diffraction gratings (photonic crystal) with the variable period one can sig- 
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nificantly increase radiation output. It is also shown that diffraction gratings 
(photonic crystal) can be used for creation of the dynamical wiggler with vari- 
able period in the system. This makes possible to develop double-cascaded 
FEL with variable parameters changing, which efficiency can be significantly 
higher that of conventional system. 

2 Lasing equations for the system with a diffraction grating (pho- 
tonic crystal) with changing parameters 

In general case the equations, which describe lasing process, follow from the 
Maxwell equations: 



here E and H are the electric and magnetic fields, j and p are the current 
and charge densities, the electromagnetic induction Di(f,t') = f £u(r,t — 
t')Ei(r,t')dt' and, therefore, Di(f,uj) = eu{r,u)Ei{r,oS), the indices i,l = 
1, 2, 3 correspond to the axes x, y, z, respectively. 

The current and charge densities are respectively defined as: 



where e is the electron charge, v a is the velocity of the particle a (a numerates 
the beam particles), 





j(r,t) = e^2v a (t)8(f- f a (t)), p(r,t) = e^5(r - f a (t)) 



(2) 



a a 
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2 i — * — * 

here j a = (1 — ^f)~ 5 is the Lorentz-factor, E(f a (t),t) (H(f a (t),t)) is the 
electric (magnetic) field in the point of location f a of the particle a. It should 
be reminded that the equation (3) can also be written as [?]: 

d-f a v a f 1 



where p Q is the particle momentum. 



Combining the equations in (1) we obtain: 



1 ()-/') \-ifi 

AB + V (V^) + ^ = -^. (5) 



The dielectric permittivity tensor can be expressed as e{r) = 1 + x{r), where 
x(r) is the dielectric susceptibility. When x <C 1 the equation (5) can be 
rewritten as: 

Ajg(f ' v-hwl 1 _ t>] ^ t>)dt> = 47r {y^r + ^ t] ) - (6) 

When the grating is ideal x(r) = ^ r Xr{f)^ %Tr ■, where r is the reciprocal lattice 
vector. 

Let the diffraction grating (photonic crystal) period is smoothly varied with 
distance, which is much greater then the diffraction grating (ptotonic crystal 
lattice) period. It is convenient in this case to present the susceptibility x(r) 
in the form, typical for theory of X-ray diffraction in crystals with lattice 
distortion [18]: 

*(0 = E ei * rW *r(r), (7) 



r 



where 3> r (f) = / f{f')dl', r(r') is the reciprocal lattice vector in the vicinity 
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of the point r' . In contrast to the theory of X-rays diffraction, in the case 
under consideration \ T depends on r. It is to the fact that x T depends on 
the volume of the lattice unit cell Q, which can be significantly varied for 
diffraction gratings (photonic crystals), as distinct from natural crystals. The 
volume of the unit cell Q(r) depends on coordinate and, for example, for a 
cubic lattice it is determined as Q(r) = , , ^ , where d; are the lattice 

periods. If x T {r) does not depend on r, the expression (7) converts to that 
usually used for X-rays in crystals with lattice distortion [18]. 

It should be reminded that for an ideal crystal without lattice distortions, the 
wave, which propagates in crystal can be presented as a superposition of the 
plane waves: 

oo 

E(r,t) = Y / Are l(k ^- U)t) , (8) 

T=0 

where k T = k + f. 

Let us use now that in the case under consideration the typical length for 
change of the lattice parameters significantly exceeds lattice period. This pro- 
vides to express the field inside the crystal with lattice distortion similarly (8), 
but with A? depending on r and t and noticeably changing at the distances 
much greater than the lattice period. 

Similarly, the wave vector should be considered as a slowly changing function 
of coordinate. 

According to the above let us find the solution of (6) in the form: 

E(r, t) = Re j £ A^iM?)-"*) j , (9) 
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where 4> T (r) = Jq k(r)dl + $ T (r), where k(r) can be found as solution of the 
dispersion equation in the vicinity of the point with the coordinate vector r, 
integration is done over the quasiclassical trajectory, which describes motion 
of the wavepacket in the crystal with lattice distortion. 

Let us consider now case when all the waves participating in the diffraction 
process lays in a plane (coupled wave diffraction, multiple-wave diffraction) 
i.e. all the reciprocal lattice vectors r lie in one plane [21,20]. Suppose the 
wave polarization vector is orthogonal to the plane of diffraction. 

Let us rewrite (9) in the form 

E(f, t) = eE(r, t) = e Re {A^W-ut) + ^(M^t) + j j ( 1Q ) 

where 

r 

0i (r) = / h(r')dl, 
o 

r r 

<hif) = j k 1 {r')dl + j r(f')dl 
o o 

Then multiplying (6) by e one can get: 

AE(r, t) - 1^ / e(r,t- t')E(r,t')dt' = Ane(^^^- + Vp(f,t)^) .(13) 

Applying the equality AE(f,t) = V(VE) and using (10) we obtain 
A^e^ 1 ^-^) = e l(Mr ^ t) [2iV(f) 1 VA 1 + iVk 1 {r)A 1 - k 2 1 (f)A 1 ], (14) 

Therefore, substitution the above to (13) gives the following system: 



(11) 
(12) 
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^e i{M " ) - UJt) [2ik 1 (r)VA 1 + iVk l {r)A 1 - k 2 1 (f)A 1 + 



u) 



2 



. 1 du 2 e (uJ, r) dAi u 2 . 1 duj 2 e^ T (uo, r) dA 2 



+ -e (u, f) Al + i- ^' — + r)A 2 + i- =^-^1 + 



i 1 <9j(r, £) \ 
+ conjugated terms = Aire | — — — h \/p(r,t) ] , 



I e *(fc(O-«0[2i£ 2 (f)VA 2 + iVk 2 (r)A 2 - k\{r)A 2 + 



2 



. 1 du 2 E (uJ, r) dA 2 ^ 2 . 1 du 2 e T (uo,r) dA 1 . 



/ i dj(r, t) \ 
+ conjugated terms = Aire I — — h Vp(r, t) j , (15) 

where the vector k 2 {r) — k\{r) +r, e (ui,r) = 1 + Xo(j)> here notation Xo(r) — 
Xr=o(^*) is used, £ r (u;,r) = Xt(t)- Note here that for numerical analysis of 
(15), if xo ^ 0, it is convenient to take the vector ki{r) in the form ki(r) = 
n\jk 2 + ^-Xo(r)- 

For better understanding let us suppose that the diffraction grating (photonic 
crystal lattice) period changes along one direction and define this direction as 
axis z. 

Thus, for one-dimensional case, when k((r)) = (k±,k z (z)) the system (15) 
converts to the following: 

\e^^'W[2ik u {z)?£ + - [k\ + k\ z (z))A x + 

uj 2 1 du} 2 eo(uj, z) dAi uj 2 1 du 2 e^ T (uj, z) dA 2l 

+—e {u, z)A l + i— — + — e- T {uj, z)A 2 + i— — + 

c 2 c 2 OLo at c 2 c 2 ouj at 

( 1 djir, t) - _ \ 

+ conjugated terms = Anel — — - 1 h Vp(f,t) , 

\c 2 at J 

l !( ^r i+fe (,)^) [2 ^ 22( ^ + MUA ^ _ 2 + k 2 {z))A2 + 

2 oz oz 

uj 2 . , A .1 du 2 e (uj, z) dA 2 uj 2 . . , .1 duj 2 e T (uj, z) dA u 



uj 2 . 1 du 2 E (oo, z) dA 2 u 2 . 

-e (u, z)A 2 + ,- Qu + -eM z)A x 

( 1 dj(r, t) \ 

+ conjugated terms = 47re — — — - h V p(r,t) , (16) 

\c 2 at J 
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Let us multiply the first equation by e -*( fc ±r ± +<^(z)-^) anc j ^ seconc [ Dy 
e -i{k±r±+4>2z(z)-uit) _ This procedure provides to neglect the conjugated terms, 
which appear fast oscillating (when averaging over the oscillation period they 
become zero). 

Considering the right part of (16) let us take into account that microscopic 
currents and densities are the sums of terms, containing delta-functions, there- 
fore, the right part can be rewritten as: 



<^r ±+Mz )-u,t) Ang f 1 + Vp(f,t)^ = (17) 

\c z dt J 
l7riUje eJ2Ut)S((r) ~ lr) a (t))e-^- + ^-^ 9(t - t a ) 6{T a - t) 



c 2 



here t a is the time of entrance of particle a to the resonator, T a is the time 
of particle leaving from the resonator, functions in (ref5) image the fact 
that for time moments preceding t a and following T a the particle a does not 
contribute in process. 

Let us suppose now that a strong magnetic field is applied for beam guiding 
though the generation area. Thus, the problem appears one-dimensional (com- 
ponents v x and v y are suppressed). Averaging the right part of (18) over the 
particle positions inside the beam, points of particle entrance to the resonator 
r±oa and time of particle entrance to the resonator t a we can obtain: 

e -*(^+*i«M-^)47re (^^|^ + Vp(r,f)) = 

= _ 4*iu>p#Mt)e l r d2f l f - i( , l( ^ A io) + Sxfu,-« t ) (ftb = 
cr S J T J 



= _ 4niup#iu(t)e << e -i {M ^ ± ,t,t ) + k ± ? ±0 -u,t) dto >>; (lg) 
c 2 

where p is the electron beam density , u(t) is the mean electron beam velocity, 
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which depends on time due to energy losses, #1 = y 1 — ^i c 2 , /9 2 = 1 — ^, 
<< >> indicates averaging over transversal coordinate of point of particle 
entrance to the resonator r_i_ 0a and time of particle entrance to the resonator t a . 

According to [22] averaging procedure in (18) can be simplified, when consider 
that random phases, appearing due to random transversal coordinate and time 
of entrance, presents in (18) as differences. Therefore, double integration over 
d 2 f± dt can be replaced by single integration [22] . 

The system (16) in this case converts to: 



uo 2 , , . .1 du 2 e (uj, z) dA 1 u 2 
. 1 du 2 e- T (uj, z) dA 2 2uj 

uj 2 . , . . 1 duj 2 e (uj, z) dA 2 uj 2 , s . 
+-s (u, z)A 2 + + -sM z)A 1 + 

. 1 du 2 s T (uj,z) dAi .2a; 

+l 7 2 du -ar = ^(M^)), 

where the currents Ji, J2 are determined by the expression 



2w 

J m = 2<Kjd m J JI__E( e -^(M,P) + e -«MM,-P)) dpj m = i ?2 (20) 




1-^- ^ = 1-^ 



2 



j = en v is the current density, A 1 = A T=0 , A 2 = A T , k± — k T=0 , k 2 = k x + r. 
The expressions for J\ for k\ independent on z was obtained in 



When more than two waves participate in diffraction process, the system (20) 
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should be supplemented with equations for waves A m , which are similar to 
those for A 1 and A 2 . 

Now we can find the equation for phase. From the expressions (11,12) it follows 
that 

d 2 (p m | 1 dv d(j} m _ dk m | k m d 2 z 
dz 2 v dz dz dz v 2 dt 2 ' 

Let us introduce new function C(z) az follows: 



C m (z)e-So'k&*' = -^C m (z), (22) 
{z>) 



dz "' ' v(z) 

z 

<f>m(z) = 0m(O) + / -^-C m (z>)dz> 

J viz 



Therefore, 



dC m (z) v(z) ( dk m k m d z 



+ ^ ■ (23) 



dz vo \ dz v 2 dt 2 

In the one-dimensional case the equation (4) can be written as: 



drz c& 

lF-^J^ ReEM ' (24) 

therefore, 

dC m (z) _ v(z) dk m t k m ed m D „ r4 +f„wj<t> m (z,t(z), P )- 



+ 7, »/ Z , , Re{A m (z,t(z))e^ z ^}, (25) 



dz vo dz vov(z) m / y 3 (z,t(z),p) 

d<f) m (t,z,p) u 

j \z=Q — K mz , (pmV'i z iP)\z=0 — Vi 

dz V 
A\\ Z =L = Ei, A 2 \ Z =L = E®, 

A m \ t=0 = 0, m = 1,2, 
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t > 0, z E [0, L], p E [— 2tt, 2tt], L is the length of the photonic crystal. 

These equations should be supplied with the equations for 7(2;, p). It is well- 
known that 

mc 2 $ = evE. (26) 
Therefore, 

d J^hA = £ ^Re{E M*, t(z))<*MM}. (27) 

The above obtained equations (20,23,25,27) provide to describe generation 
process in FEL with varied parameters of diffraction grating (photonic crys- 
tal). Analysis of the system (25) can be simplified by replacement of the 
j(z,t(z),p) with its averaged by the initial phase value 



Note that the law of parameters change can be both smooth and stair-step. 
Use of photonic crystals provide to develop different VFEL arrangements (see 



It should be noted that, for example, in the FEL (TWT,BWO) resonator 
with changing in space parameters of grating (photonic crystal) the electro- 
magnetic wave with depending on z spatial period is formed. This means that 
the dynamical undulator with depending on z period appears along the whole 
resonator length i. e. tapering dynamical wiggler becomes settled. It is well 
known that tapering wiggler can significantly increase efficiency of the undu- 
lator FEL. The dynamical wiggler with varied period, which is proposed, can 








Fig.l). 
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Fig. 1. An example of photonic crystal with the thread arrangement providing 
multi-wave volume distributed feedback. Threads are arranged to couple several 
waves (three, four, six and so on), which appear due to diffraction in such a struc- 
ture, in both the vertical and horizontal planes. The electronic beam takes the whole 
volume of photonic crystal. 

be used for development of double-cascaded FEL with parameters changing in 
space. The efficiency of such system can be significantly higher that of conven- 
tional system. Moreover, the period of dynamical wiggler can be done much 
shorter than that available for wigglers using static magnetic fields. It should 
be also noted that, due to dependence of the phase velocity of the electromag- 
netic wave on time, compression of the radiation pulse is possible in such a 
system. 

3 Conclusion 

The equations providing to describe generation process in FEL with varied 
parameters of diffraction grating (photonic crystal) are obtained. It is shown 
that applying diffraction gratings (photonic crystal) with the variable period 
one can significantly increase radiation output. It is mentioned that diffrac- 
tion gratings (photonic crystal) can be used for creation of the dynamical 
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wiggler with variable period in the system. This makes possible to develop 
double-cascaded FEL with variable parameters changing, which efficiency can 
be significantly higher that of conventional system. 
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